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Abstract 

Results are presented on the discovery potential for MSSM neutral Higgs bosons 
in the mh — max scenario. The region of large tan/?, between 15 and 50, and mass 
between 95 and 130 GeV is considered in the framework of the ATLAS experiment 
at the Large Hadron Collider (LHC), for a centre-of-mass energy ^/s = 14 TeV. This 
parameter region is not fully covered by the present data either from LEP or from 
Tevatron. 

The h/A bosons, supposed to be very close in mass in that region, are studied in 
the channel h/A — > fi~^fi~ accompanied by two b-jets. The study includes a method 
to control the most copious background, Z — * accompanied by two b-jets. A 

possible contribution of the H boson to the signal is also considered. 



1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) is the most investigated extension of 
the Standard Model (SM). 

The theory requires two Higgs doublets giving origin to five Higgs bosons: two CP-even 
neutral scalars, h and H (h is the lighter of the two), one CP-odd neutral scalar, A, and one 
pair of charged Higgs bosons, [1^3]. The discovery of any one of these particles is a crucial 
element for the confirmation of the model. This is a key point in the physics program of future 
accelerators and in particular of the LHC. 

After the conclusion of the LEP program in the year 2000, the experimental limit on the 
mass of the Standard Model Higgs boson H was estabhshed at 114.4 GeV with 95% CL [4]. 
Limits were also set on the mass of neutral [5] and charged [6] MSSM Higgs bosons for most of 
the representative sets of model parameters. 

The motivation of this study is to explore the potential of the ATLAS detector for the 
discovery of neutral MSSM Higgs bosons in the parameter region not excluded by the LEP and 
Tevatron data. We shall focus on the search for h, the hghtest of the neutral Higgs bosons. Its 
mass, taking account of radiative corrections, is predicted to be smaller than 140 GeV, see [5] 
and references therein. The search for a mass close to the mass of the Z boson, mz, will be a 
challenging test of detector performance and of the analysis method of disentangling the signal 
from the background. 

In the first part of this paper we review the MSSM framework, the production mechanism 
in hadron collisions, the present experimental situation and the discovery potential at the LHC. 

In the second part, after describing the Monte Carlo generator and the software tools used, 
we discuss the detector performance relevant for this search of h, A and H, the analysis strategy 
and the results of the scan over the MSSM {ttia, tan/3) plane. Details of this analysis are given 
in Ref. [7] and Ref. [8]. 

In the conclusion, results are presented on the neutral MSSM Higgs bosons discovery po- 
tential at the LHC based on the ATLAS detector. 

2 Minimal Supersymmetric Standard Model 

We discuss a few points of the model, useful for the present analysis. For a complete review 
see Refs. [2,9]. At tree level, the masses of the five Higgs bosons of the MSSM are related by 
the following equations: 



where mw and mz are the W and Z masses, respectively, and tan/3 is the ratio of the vacuum 
expectation values of the two Higgs fields. 

The MSSM model may be constrained by the assumption that the sfermions (scalar fermions) 
masses, the gaugino masses and the trilinear Higgs-fcrmion couplings must unify at the Grand 
Unification scale (GUT). In one of the possible constrained models the parameters chosen are: 

• MsusY, a common mass for all sfermions at the electroweak scale. 

• M2, a common SU{2)l gaugino mass at the electroweak scale. 
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• /i, the strength of the supersymmetric Higgs mixing. 

• tan/?, the ratio of the vacuum expectation values of the two Higgs fields . 

• A = At= Ab a common trilinear Higgs-squarks coupling at the electroweak scale. It is 
assumed to be the same for up-type squarks and for down-type squarks. 

• m-A, the mass of the CP-odd Higgs boson. 

• mg, the gluino mass. 

Three of these parameters define the stop and sbottom mixing parameters Xt = At— ^ cot/3 
and Xb = Ab — /i cot/3. 

Whereas the particle spectrum depends on all the parameters mentioned above, the Higgs 
sector depends, at tree level, on only two parameters which can be taken to be tan/3 and mA, 
as in Eq. ([1]). The other parameters only enter through radiative corrections, but change the 
mass prediction of Eq. ([T]) (where it is limited to m-h < mz), allowing the mass of h to reach 
higher values (~ 130 GeV in some scenarios, see Sec. 14.11) . 

Among all possible CP-conserving benchmark scenarios the so-called mh — max scenario 
(Table [1]) has been considered [5]. This scenario corresponds to the maximum value of the stop 
mixing parameter Xt=A - /i cot/5 = 2Msusy- Here the theoretical bound on the mass of the 
h is highest (hence the scenario's name) and experimental limits are less constraining. Also, 
the range of excluded tan/5 values for given values of rrit (the top mass) and Msusy is the most 
conservative one. 



Parameter 


mh — max 


Msusy [GeV] 


1000 


/i [GeV] 


-200 


Ms [GeV] 


200 


Xt=A- fx cot f3 


2MsusY 


mg[GeV] 


O.SMsusY 


mA[GeV] 


0.1-1000 


tan/3 


0.4-50 



Table 1: CP-conserving benchmark mh — max scenario . 



We shall focus on the mh — max scenario as the most promising for the search of the h 
boson, referring to it hereafter as MSSM. 

At tree level the MSSM Higgs boson couplings to fermions and massive gauge bosons are 
obtained from the SM Higgs boson couplings via correction factors [10]. They depend on the 
parameters (3 (already introduced) and a, the mixing angle which diagonalizes the CP-even 
Higgs boson mass matrix. The two parameters are related by the following expression: 

^ m^ — m^ 

cos 2a = - cos 2(3 -4 1 (2) 

mj^-mi 

At high tan/? the MSSM correction factors to the SM Higgs bosons couplings to fermions 
and massive gauge bosons are larger for down-type quarks (b) and leptons (r and /i) than for 
up type-quarks. Thus, the associated bbh production is enhanced and becomes the dominant 
process in the production of h bosons in the high tan/5 region. 
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The Feynman diagrams contributing to the process gg^ bbh hhn^n and qq-^ bbh 
bb/i+/i~ are shown in Fig. [H 




Figure 1: Typical diagrams contributing at "tree level" to the process gg^ bbh bb/i+/i and qq^ bbh 
bb/x+yit^. 

We shall consider mainly the decays to jX^jJi~ . Indeed, although the Higgs boson couplings 
are proportional to the fermion mass, thus resulting in a branching ratio to r+r^ higher than 
to by a factor (^)^, the experimental conditions favor the fi'^fi~ channel^. 

In the region of high tan/3 and mh ~ 100 GeV the CP-odd neutral Higgs boson A has a mass 
only slightly higher than the CP-even h and a competitive branching ratio for the fi~^fi~ decay 
channel [10]. Also, cross-sections and widths, which are functions of the parameters tan/3 and 
TJiA, are close in some regions of the parameter space (Sec. 14.11 ). Thus, in these regions the h 
and A bosons are indistinguishable from an experimental point of view, and it is more correct 
to think in terms of a h/A search. In the following we refer to the boson sought as the h/A 
boson (however its mass is noted mh or mA, accordingly) . The degeneracy between h and A 
is less pronounced near the higher mass limit of h . 

3 Experimental search for Minimal Supersymmetric Stan- 
dard Model Higgs 

3.1 LEP and Tevatron results 

High precision tests of the Standard Model have been performed at LEP setting a combined 
limit of mH> 114.4 GeV for the mass of the SM Higgs boson [4]. 

Again at LEP, the validity of the Minimal Supersymmetric Standard Model has been inves- 
tigated within the constrained framework of Sec. [2l For the mass of the charged MSSM Higgs 
bosons a combined limit mH±> 78.6 GeV was obtained [6]. Searching for neutral CP-even and 
CP-odd MSSM Higgs bosons, no indication of signal was found up to a center-of-mass energy 

The production advantage of the r+r" channel is counterbalanced by the difficulty of identifying the 
hadronic decay of a T-jet in hadronic events, by a smaller acceptance of the detector and by a worse mass 
resolution due to the presence of neutrinos in the final state. Instead, with a final state like h A*^/^^ ATLAS 
would exploit the excellent combined performance of the muon spectrometer and inner detector. 
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(GeV/c^) nij^ (GeV/c^) 

Figure 2: The combined LEP results for the search for the MSSM neutral Higgs bosons (from Ref. [5]). The 
figure shows the theoretically inaccessible regions (light-grey/yellow) and the regions experimentally excluded by 
LEP searches, at 95% C.L. (medium-grey/light-green) and 99.7% C.L. (dark-grey/dark-green), for the nih-max 
scenario with the top mass mt= 174.3 GeV, in two projections of the MSSM parameters (mh, "T'a), (m^h, tan/3). 
The dashed fines indicate the boundaries of the regions which are expected to be excluded, at 95% C.L., on 
the basis of Monte Carlo simulations with no signal. In the (mh, tan/3) projection, the upper boundary of the 
parameter space is indicated for four values of the top mass; from left to right: mt= 169.3, 174.3, 179.3 and 
183.0 GeV. 



of 209 GeV [5]. The corresponding lower limits on the masses were set as a function of tan/3 for 
several scenarios. In the mi^-'max scenario (Fig. [2]) with a top mass mt=174.3 GeV the limits 
for tan/? > 10 at 95% CL are approximately : 

'Tih, > 93 GeV 

A complementary search, providing sensitivity in the region tan/5 > 50 has been performed 
at the Tevatron Collider at ■\/i= 1.96 TeV. In the MSSM scenario, a significant portion of the 
parameter space has been excluded by the DO Collaboration, down to tan/5 = 50 as a function 
of TTT-A, by studying the associated production with two b quarks of h/A/H bosons and their 
decay into bb [11]. Comparable results have been obtained by the CDF Collaboration exploring 
the h/A/H decays to r+r^, but extending the excluded region to higher values of niA [12]. 

3.2 LHC discovery perspectives 

The LEP and Tevatron data don't exclude the parameter space defined by tan/3 larger than 
10 and smaller than 50. Therefore, a natural continuation of the LEP and Tevatron physics is 
the investigation of the possible existence of MSSM Higgs bosons in this region of tan/?. The 
ATLAS [13] and CMS [14] experiments starting in the near future at the Large Hadron Collider 
(LHC), at CERN, constitute an excellent laboratory for such search. 

The prospect for the detection of MSSM Higgs bosons at LHC was evaluated for benchmark 
sets preventing Higgs boson decays to SUSY particles [10, 13] and focusing on the discovery 
potential of decay modes common to MSSM and SM Higgs bosons [13]. It was concluded that 
the complete region of parameter space tjia = 50 - 500 GeV and tan/5 = 1 - 50 is open to Higgs 
boson discovery by the ATLAS experiment, already with an integrated luminosity of j Cdt = 
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30 fb~^, and that over a large part of this region more than one Higgs boson and more than 
one decay mode could be observed - the detection of a signal in more than one decay channel 
would constitute strong evidence for the MSSM model. It was also found that the region in the 
(mA, tan/5) plane which corresponds to rrih ^ 100 GeV and tan/5 > 10 is only accessible by a 
neutral h/A boson decaying to yU'^/i" or r+r^ [10, 13], and by a charged boson decaying to 
riy [15]. 

More recently the h boson discovery potential in the MSSM scenario has been investigated 
[16] at two luminosities, J C dt= 30 fb^^ and J C dt= 300 fb^^ . At low luminosity the h — * t~^t~ 
decay mode represents the main contribution to the discovery potential and covers most of the 
parameter space not yet explored. However the contribution of bbh — > /i"*"/!" appears to be 
crucial in the region of moderate tan/3 and mass close to mz- The channel bbh — /i"*"//" which 
requires an excellent performance in fi detection and b-tagging is well suited to the ATLAS 
experiment thanks to a design giving high performance from the muon spectrometer and the 
inner detector. 

At high luminosity channels such as: h 77, h ^ ZZ — > 4£ and h — bb in associated 
production with tt give a significant contribution. The channel h —>■ 77, which requires an 
excellent M^^ mass resolution and jet/7 separation, corresponds to MSSM rates suppressed 
with respect to the SM case but for a limited region of the parameter space where they could 
even be slightly enhanced. As for the channel h bb, only the tth production followed by the 
h — > bb decay can be observed clearly above the background, thus the extraction of the signal 
requires the identification of four b-jets and an excellent b-tagging performance. In the MSSM 
case the rates could be enhanced by 10-20% over the SM rates. 

3.3 Background processes 

The prospects for the detection of MSSM Higgs bosons at LHC depend heavily on the suppres- 
sion of the background sources: 

• Z/7* production with two b-jets and a subsequent decay into a pair. The cross 
section for this process is a few orders of magnitude larger than for the signal. As an 
example, we quote for the Z boson a^^^z ■ Brz^^+^- ^ 22.8 pb and for each of the 
h and A bosons cTh/Abb ' ^^h/A^/i+ti- ~ 0.24 pb (at tan/5 = 45 and rrih = 110 GeVjf^. 
The corresponding diagrams are shown in Fig. [3] and Fig. [T|, respectively. It is clear 
that they differ only in the kind of boson produced, h (or A) in the signal and Z/7* 
in the background. The distinction between the signal and the background when mh is 
approaching mz becomes then an extremely hard task due to the similar topology of the 
decays [7], although angular distributions differ somewhat due to the fact the boson is a 
scalar in one case and a vector in the other. 

• Z production with two jets not originating from b-quarks. The cross section is ~ 24 
times that for Z with two b-jet£I] and can contribute to the background in case of jet 
misidentification. An estimate of a possible impact on the significance of this analysis is 
reported in Sec. 17.21 

• ZZ associated production, when one Z decays into bb and the second one decays into 
/i^/i". This process has a cross section of the same order of magnitude as the signal, 

2) Evaluated from AcerMC(2.3) [17] and PYTHIA 6.226 [18] (with mz > 60 GeV). 

3) Evaluated from PYTHIA 6.226 [18]. 

'^^Evaluated (for this purpose) from bbZ and jjZ cross sections by SHERPA 1.0.9 [19] (with mz > 60 GeV). 
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Figure 3: Typical diagrams contributing at "tree level" to the process gg^ Z/7*bb bb and 

qq — > Z/7*bb — s- /x+/^~bb. 

(Tzz • Brz^^+^- • Brz^bb ~ 0-13 pb but can be easily suppressed using the kinematic 
characteristics of the events, see following sections. 

• tt associated production followed by a top-quark decay into a b-quark and a W boson and 
a subsequent W decay in fiu. The cross section of this process is ati ■ Brt^bw ■ Brw^/^i/ • 
Brt^bw ■ Brw^/^j/ ~ 5.71 pb[^. The presence of two neutrinos implies missing transverse 
energy in the event (see following sections) thus allowing this background to be strongly 
reduced. One would want to discriminate the signal from the background on the basis of 
the different b-jets characteristics (the two background b-jets are usually more energetic 
than those accompanying the signal and the probability of their identification is higher), 
but the requirement that the two b-jets be identified will suppress the signal more than 
the background. 

4 Monte Carlo simulation 

We improve on the analyses reported in Sec. 13.21 with a full Monte Carlo simulation of the 
experiment (data generation, reconstruction and analysis). The exploration of the unexcluded 
MSSM parameter space, with a view to either discovering a supersymmetric Higgs boson or 
excluding the model considered, constitutes the motivation of the analysis described in this 
paper. 

The considerations of Sec. [2] and Sec. [3] suggested a search for h/A — »■ decays accom- 

panied by two b-jets. The extraction of the h/A boson signal from the competing enormous 
background of Z decays in the mass region close to constitutes a challenging search, where 
all performances of the experimental setup have to be exploited. 

To this purpose we have generated signal and background Monte Carlo events using the 
PYTHIA program (v. 6. 226) [18] as specified below, for a center-of-mass energy ^/s= 14 TeV. 
The efficiency of the selection criteria, the detector acceptance and the purity of the data sample 
are estimated from these events and the ATLAS detector response [13]. The latter is simulated 
using the GEANT program [20,21] which takes into account the effects of energy loss, multiple 

5) Evaluated from PYTHIA 6.226 [18]. 

6) Evaluated from PYTHIA 6.226 [18]. 
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scattering and showering in the detector through an interface called ATHENA (v. 10. 0.1). 

The number of events used in this analysis corresponds to an integrated luminosity j Cdt = 
300 fb~^ with the exception of the channel bbZ hhfi~^fi~. In this latter case, for practical 
reasons, a number of events corresponding to half the mentioned luminosity has been generated. 
We note that JCdt = 300 fb~^ corresponds to the integrated luminosity expected after three 
years of data taking. 

4.1 Signal 

The MSSM neutral Higgs bosons h, A and H were generated in associated production with two 
b-quarks using the PYTHIA program (v.6.226) through the ATHENA interface (v. 9. 0.4). The 
parameters of the model were given the values [22] shown in Table [2] (for tan/3 and the 
range of values scanned is shown). In this parameter region, as mentioned in Sec. Ej the cross 
section, mass and width of h and A bosons are close while the H cross section is one order of 
magnitude lower than the h/A cross section. 

As is done in the Pythia code, we have taken ttia as input parameter. The values of mh 
and mn are then derived from PYTHIA as a function of rrix- 



Parameter 




value 


common gaugino mass 


M2= mi/2 [GeV] 


200 


gluino mass 


mg [GeV] 


800 


strength of supersymmetric Higgs 


/i [GeV] 


- 200 


ratio of Higgs fields 


tan/5 


15-50 


common scalar mass 


mo [GeV] 


1000 


squark left 3 gen 


[GeV] 


1000 


sbottom mass 


^bR [GeV] 


1000 


stop mass 


[GeV] 


1000 


stop-trilinear coupling 


A 


2440 


mass of CP-odd boson 


niA [GeV] 


95-135 



Table 2: Parameters for A, h and H generation in PYTHIA (v.6.226) [18]. 

The signal bbh — >■ hhn~^fi~ has been simulated in 104 points of the parameter space (mA, 
tan/?), corresponding to eight steps in tan/?, chosen equally spaced between 15 and 50 and 
thirteen steps of 2.5 GeV in mh between 95 GeV and 125 GeV (the largest value allowed 
by PYTHIA (v.6.226)). These are also the points where the decay bbA bb/i^/i~ has been 
simulated. 

The values of the mass, cross section and width of the A/h and H neutral bosons are reported 
in Ref. [7] for all the points analyzed of the (mA, tan/?) plane. For the following discussion a 
reference point has been chosen at tan/? = 45, mA=110.31 GeV (mh= 110 GeV, mH= 127.46 
GeV). The corresponding width is rA= 4.28 GeV (rh= 4.20 GeV, 0.05 GeV), while the 
production cross section times the branching ratio for the decay to is o"bbA^bb/i+At- = 0.243 
Pb (t^bbh^bbM+M- = 0-245 pb, abbA^bbM+M- = 0-0016 pb). 
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4.2 Background 

For convenience we list again below the sources of background considered in Sec. 3.3 which 
were fully simulated. The Z production associated with two jets not originating from b as 
yet has not been implemented in the AcerMC2.3 code and has been ignored at this level (see 
however end of Sec. 17.21) . 

• bbZ —>■ bbyU+yU". The events are generated by AcerMC(2.3) [17] with the hadronization 
process of PYTHIA (v.6.226) 111 At the generator level a cut-off is applied to the 
invariant mass, M^J^J > 60 GevB. 

• tt ^ hh^^fi^uu. The events are fully generated with PYTHIA (v.6.226). 

• ZZ ^ bbyU^/i~. The events are fully generated with PYTHIA (v.6.226). 



process 
(background) 


[Pb] 


]\^exp3oo 




w 


bbZ bb/i+yU~ 


22.789 


6836700 


3314000 


2.06 


tt bbyU^/i~z/P 


5.71 


1713420 


1806437 


0.95 


ZZ —> hhfi~^fi~ 


0.1273 


33819 


97244 


0.35 



Table 3: Background cross section times branching ratios, Cbb/i+/j-j number of expected events for JCdt = 
300 fb"\ N'^'^P^o, number of Monte Carlo generated events, N^^, and their weight in the analysis, w, for the 
three processes considered. 

The cross section times branching ratios, together with the number of expected events for 
JCdt = 300 fb~^, the number of Monte Carlo generated events and their weight, is reported 
in Table [S] for the three background processes. 

Among these the main contribution comes from bbZ hh^'^ fi~ , and is affected by a cross 
section uncertainty arising from uncertainties on QCD, QED couplings (~ 10% [17]) and from 
higher order corrections (~ 25% [23]). As we intend to account for this large uncertainty with a 
data-driven method (Sec. 16. 3|) . a sample of bbZ —>■ e"^e~ events, ~ 600000 events, corresponding 
to an integrated luminosity JCdt = 30 fb~^ with a cross section a ~ 22.8 pb, has also been 
simulated. Owing to the large number of events simulated, there is only a minor statistical 
uncertainty associated with the evaluation of the tails of Z — > /x^/x" background. 

5 The ATLAS detector 

Here we shall summarize the main features of the ATLAS detector [24] which are relevant for 
the present analysis. For a more detailed information we refer to [25] : 

• Magnet System: This consists of a solenoid providing a 2 T-m bending power in the 
inner detector, and a barrel air toroid completed by two end-cap toroids, with a typical 
bending power of 6 T-m and 3 T-m respectively. 

^)Using PHOTOS package for inner bremsstrahlung generation. 

A low energy cut is fixed on inner bremsstrahlung photons at Pt > 5 GeV. 
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• Inner Tracking System: It has been designed to measure as precisely as possible 
and with high efficiency the charged particles emerging from primary interactions in a 
pseudorapidity range |?7| < 2.5. It is installed inside in the solenoid magnetic field and 
consists of Si pixels and silicon strip detectors near the interaction point, and strawtubes. 
The expected transverse momentum resolution of the inner tracker for a 100 GeV charged 
particle at |?7| = is ^ 3.8%. 

• Electromagnetic Calorimeter System: The excellent energy resolution and particle 
identification for electron, photons and jets demanded from physics is realized from a 
liquid argon-lead sampling calorimeter with accordion shape in the barrel and end-cap 
regions. The energy resolution expected is ^ = ^^'J^'^" © 0.35. 

• Hadronic Calorimeter: It is a copper-liquid argon calorimeter in the end-cap region 
and a Fe-scintillator calorimeter scintillators in barrel region. The liquid argon tungsten 
forward calorimeters extend the coverage to |?7| =4.9. The expected hadronic jet energy 
resolution is | = ^ © ^ © 0.013. 

• Muon spectrometer. The reconstruction of muons at highest luminosity is one of the 
most important point in ATLAS design. The ATLAS toroidal magnet field provides a 
muon momentum resolution that is independent of pseudorapidity. The spectrometer 
is constitued by: a) The precision tracking chambers made of Monitored Drift Tubes 
(MDT) covering the rapidity range |?7| < 2.7, b) The Cathode Strip Chambers (CSC) and 
transverse coordinate strips cover the most forward rapidity region ( Ir^l = 2.0 — 2.7) in 
the inner most stations of the muon system, c) Resistive Plate Chambers (RPC) in the 
barrel region ( |?7| < 1) and Thin Gap Chambers (TGC) in the end-cap region provide 
muon triggers and measure the second coordinate of the muon tracks. The expected 
momentum resolution ranges from about 1.4% for 10 GeV muons to 2.6% for 100 GeV 
muons atlr^l =0 

6 Preliminaries to Monte Carlo data analysis 

Two points are crucial for our analysis: 

• The muon reconstruction efficiency in the analysis acceptance and the invariant 
mass resolution. 

• The b-jet identification. 

6.1 Efficiency and resolution studies with bbZ — > bb/z"'"/x~ events 

The /i reconstruction performance of the apparatus was studied with a sample of bbZ — >■ 
bb/i"*"/!" events [26]. For this purpose only events with two reconstructed muons of opposite 
charge found within the \ri\ < 2.5 acceptance were considered. The contribution of tracks 
mimicking muons at reconstruction level was found to be negligible in a Z^fi'^fi" sample [27]. 
The pileup effect seems to be also negligible based on the information presently available [28]. 

The distributions of the transverse momentum, Pt^j, pseudorapidity, 77^, and polar angle, 
(p^, of reconstructed muons (see Fig. HI green histograms) reproduce with good efficiency the 
generated data (see Fig. HI brown histograms). 
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Figure 4: Distributions of the transverse momentum Pr^,, pseudorapidity rj^, and polar angle cj)^ for muons from 
bbZ(-^ M^A* ) Monte Carlo events, generated (brown) and both reconstructed (green) within the \r]\ < 2.5 
acceptance. 



The distribution of the reconstructed dimuon invariant mass is shown in Fig. [3] (top) together 
with the result of a Gaussian fit. The sohd hne corresponds to the following values for the fit 
parameters: 

< m;";; >= 90.47 ± 0.05 GeV (r{M'^p = 3.02 ± 0.06 GeV (3) 

The fit mean is smaller than the nominal value of mz [29]. The Z natural width contributes 
to cr(M^^^) for approximately 1.9 GeV thus implying a measurement accuracy a^-es = 2.3 GeV. 
This is just the value obtained unfolding the reconstructed distribution with the distribution of 
the Z generated mass, Mgg^. When the difference between reconstructed and generated mass 
M^^^ — Mggj^ is plotted. Fig. [5] (bottom), the Gaussian fit of the distribution gives the following 
results for the fit parameters: 

< M'^^ - M^^^ >= -0.82 ± 0.03 GeV a^es = 2.35 ± 0.03 GeV (4) 

From this study we conclude that in the mz region the reconstructed fi^fj^" invariant mass 
distribution shows a mean value shifted by 820 MeV with respect to the nominal value, towards 
the low mass region, and a resolution of ~ 2.6 %. These results will be used later in Sec. 17. 2[ 

6.2 b-tagging studies 

The two b-tagging algorithms called 3D and SV2, which are in use within the ATLAS Collab- 
oration [30], were studied with a subsample of bbh bb/i"^/i~ events. 

Starting from the impact parameters (transverse and longitudinal) and their significances, 
both algorithms assign a weight to each track of the jet, that is to each track within a cone of 
opening angle AR = \/ Arj'^ + A0^, with at least one track reconstructed in the tracker. The 
track weight is the ratio between the likelihood functions for a b-jet track and for a track of a 
light jet (from a light quark, u, d, s, c or a gluon). In turn the jet weight Wjet is defined as the 
sum of the logarithms of the tracks weights. This allows b-jets (wjet > Wcut) to be discriminated 
against light jets (wjet < Wcut) - the value Wcut is chosen, depending on the physics under study. 
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Figure 5: (top) Distribution of the invariant mass -M^"^ reconstructed in bbZ(^ ) events, (bottom) 

Distribution of the difference between reconstructed and generated mass M^'™ — M^^,^ for the same sample of 
events. 
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such as to have a high efficiency eb for b-jet identification and a high rejection factor Rj for 
hght jets. 

The SV2 algorithm improves on the 3D algorithm by using additional information on other 
variables not too strongly correlated with the impact parameters (as e.g. fraction of the jet en- 
ergy and invariant mass of all particles at the secondary vertices, number of two-track secondary 
vertices) . 

In Table IH the results of a study of a bbh bb/i^/i~ sample are reported. Similar results 
have been obtained from a bbZ bb/i^;U~ sample. For both samples, a selection cut, -Prjet > 
15 GeV, is applied, since below this value the efficiency of b-identification drops. The study is 
limited to the inner detector acceptance |?7| < 2.5 and to two values of the jet opening angle, 
AR = 0.7 and AR = 0.4 



Algorithm 


AR 


Cb 


Rj 


Wcut 


AR 


Cb 


Rj 


Wcut 


SV2 


0.7 


49% 


71 


1 


0.4 


50% 


58 


1 


3D 


0.7 


55% 


27 


1 


0.4 


54% 


26 


1 


SV2 


0.7 


46% 


219 


2 


0.4 


46% 


200 


2 


3D 


0.7 


49% 


52 


2 


0.4 


49% 


50 


2 



Table 4: b-tag performance. The table shows, for the two algorithms used, the jet opening angle AR, the 
efficiency on b-jet identification eb, the rejection of light jets Rj and the value of Wcut (see text). 

As results of this study, it is possible to conclude that: 

• The SV2 algorithm, which uses more information on b decays products, provides a higher 
rejection rate with little loss in efficiency. 

• No relevant improvement is obtained by shrinking the jet cone opening angle, AR, from 
0.7 to 0.4. Therefore the ATHENA (v. 10.0.1) defauh value, AR=0.7, will be used in the 
following analysis. 

• The cut value Wcut = 1 turns out to be a good compromise with a b-jet identification 
efficiency of ~ 50% and a rejection of light jets by a factor of ~ 70. 

In case of high luminosity, possible track pileup doesn't have a significant influence on b-tagging 
efficiency for this class of events. Their topology with two muons emerging from the primary 
vertex facilitates the selection of the correct primary vertex out of many interaction vertices. 

6.3 Method for bbZ — > bb/x+/x~ background subtraction 

The method proposed in [8] exploits the two following points (at the level of particle generation): 

a) the rate of h/A e~^e~ is expected to be suppressed with respect to the signal h/A ^ f^'^f^" 

by a factor j , 

b) the rate of the background bbZ — bb/i+/i~ is equal to the rate of bbZ bbe"'"e~ because 
of the production diagrams which are the same, and of the lepton coupling universahty in the 
Z decay. 

In this context the associated Z production and decay in the channel bbZ — > fi'^fi~ has been 
studied using a control sample of bbZ e^e~ events. The effect of inner hremsstrahlung (IB) 
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radiation has been investigated and corrected for; the impact in the event reconstruction is not 
large. The ratio of the number of reconstructed events from the two samples in the region of 
mass higher than mz, interesting for new physics search, is stable and implies correction factors 
close to one [8] . 

As a result, barring corrections for different inner hremsstrahlung and detector response, 
the number of bbZ — > bbe"'"e~ gives directly the number of background events bbZ — > bb/i"*"/!". 
Details of the method are given in Appendix |Al 

6.4 Significance of the search 

Statistical methods used in Higgs boson searches are discussed in Ref. [31]. Here the significance 
of a search is given using ^ as a statistical estimator, where S is the number of signal events 
(h/A or h/A/H), and B the number of background events. Discovery means that the signal is 
larger than 5 times the background statistical error > 5). The probability of a background 

fluctuation of this size is less than ^ 2.87 ■ 10~^. A search resulting in > 3 is interpreted 
as an indication of new physics. 

7 Monte Carlo data analysis 

7.1 bbA bb^t+/x~, bbh bb/x+^t~, bbH bb/Lt+/x~ 

The signature of the h/A channel is a pair of well isolated high-energy muons with opposite 
charge and two hadronic jets containing b quarks. The invariant mass of the reconstructed 
muons is supposed to originate from a h or A boson and must be compatible, within the mass 
resolution, with the corresponding mass, mh or mp^. 

No simulation of the trigger is implemented; however two muons of transverse momentum 
above trigger experimental threshold are required and the detector acceptance considered is 
|?7| < 2.5. The event selection is divided in three steps: preselection, tt veto selection, final 
selection. 

The preselection, cuts 1-2-3 in Table [5l requires at least two well identified opposite charge 
muons with Pt^j > 10 GeV in the pseudo-rapidity range | r/^l < 2.5 (cut 1), thus satisfying 
the trigger requirement. The presence of a jet pair with -Prjet > 10 GeV and | ?7jet| < 2.5 is as 
well demanded, without any b-identification requirement (cut 2). 

A further requirement (cut 3) is that at least one of these jets be identified as originating from a 
b-quark with PTjet > 15 GeV, that is the energy lower limit for a reliable b-jet identification [30] 
(Sec. 16.21) . These cuts are designed to select events fulfilling the minimum conditions to be 
analyzed later. The events excluded are in any case not suitable to undergo any further analysis. 

The tt veto selection is designed to suppress this background, characterized by a large 
missing transverse energy due to the presence of neutrinos. E^^^^ is required to be less than 
45 GeV (cut 4). Other cuts (cut 5-6) are applied on the first and second most large muon 
transverse momentum, Pt^ii and Pt^2, and (cut 7-8) on the first and second most large jet 
transverse momentum, Prjeti and PTjet2 • At least one of these jets was previously identified as 
originating from a b-quark (cut 3). 

The final selection criteria are designed mainly to disentangle the signal events from the 
irreducible Z background. The main requirement is that the fj^~^fi~ invariant mass has to lie 
inside a window around h/A mass, determined by the natural width of the bosons and by the 
experimental resolution of the fi~^fi~ invariant mass. Sec. [61 
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The muons originating from b decays in bbbb events, with a cross section of ~ 500 pb [32], 
can mimic the signal events. To avoid this possible contamination a low hadronic activity near 
both muons is required. The isolation criteria demands that the sum of the charged track 
momenta in a cone (AR < 0.2) around the muon direction be less than 5 GeV. These selection 
cuts are summarized in Table 



Cut Number 


Variable 


1 


> 2, Ptm > 10 GeV, \r]^\ < 2.5 


2 


iVjct > 2, PTjet > 10 GeV, Ir^jetl < 2.5 


3 


A^'bjct > 1,-PTjet > 15 GeV, Wbtag >1 


4 


Ef""" < 45 GeV 


5 


25 GeV < Ptmi < 95 GeV 


6 


20 GeV < Pt/.2 < 60 GeV 


7 


PTjet2 < 40 GeV 


8 


Pxjeti < 70 GeV 


9 


IM'J'J - tbaI or iMj^J - mhl < 2a 


10 


S Prtracks < 5 GeV in cone AR < 0.2 



Table 5: Selection criteria for h/A — > fi^ fi with two b-quarks final state. 

A key point of the selection is the determination of the mass window used (cut 9). Its 
value is determined as a function of Fh (Ta), the total width of the h (A) Higgs boson, and the 
experimental mass resolution = 2.6% (Sec. 16. ip . The mass window is centered on the h (A) 
mass corrected for the bias k in the muon reconstruction, m^^^ = mh,A ~ k with k= —820 MeV 
(Sec. 16.11) . Its limits, ^j^^, are defined by: 

where f is the standard deviation factor corresponding to a chosen probability. In our analysis 
f = 2, which corresponds to including 97.7% of the signal. The event is selected if the 
invariant mass is inside the window limits either for the A or the h mass. 

For the bbH bb/i^/i" channel the procedure is the same as for the h/A search, but for the 
selection window. The window limits are obtained from Eq. [5] with mn (the H mass) replacing 
mh,A in mh°A. 

7.2 Reference point 

The analysis procedure is first discussed for the reference point m-A = 110.31 GeV, tan/5 = 45 
(mh = 110.00 GeV, mu = 127.46 GeV). For a given MSSM neutral Higgs boson (h. A, H) at 
this point. Table [6] shows the production cross section times the /x"*"/!" decay branching ratio, 
(Thhti+ii-y the expected number of signal events at fCdt = 300 fb~^, N'^'^p^oo^ number of 
Monte Carlo generated events N'^^ and their weight w. The weight of the signal events is close 
to one as it is for all (mA, tan/3) points [7]. The background weights (Table ED are also close to 
one. 

We shall justify the cuts by showing examples of distributions obtained for the bbh —>■ hhfi~^fi~ 
signal and the bbZ hh fi^fi~ background, and refer otherwise to [7] either for variables 
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different from tlie ones displayed or for a different signal (bbA —>■ hhfi^fi^ or bbH hhfi^fi^) 
or for a different background (tt —>■ hhfi^^^uu or ZZ —>■ bb/i+yU^). In all cases the reconstructed 
distributions without cuts compare to the generated distributions as expected. 



process 
(signal) 


(^hh fi+fi- 
[pb] 


jv^expaoo 




w 


bbh bbyU^/i~ 


0.245 


73500 


76517 


0.96 


bbA — > hhfi^fi^ 


0.2433 


72900 


74965 


0.97 


bbH — i> bb/i+/i~ 


0.001619 


486 


600 


0.88 



Table 6: h, A, H production at the reference point tan/3 = 45, toa = 110.31 GeV(mh = 110.00 GeV, mn = 
127.46 GeV): process, cross section times decay branching ratio, crbb^^+^j-i number of events expected for 

JCdt = 300 fb"\ ]\[<:''P3oo^ number of Monte Carlo generated events, N^'^, their weight w. 

The distributions of the number of muons A*"^ and of the muon transverse momentum Pt/x 
in the processes bbh — >• hhfi~^fi~ and bbZ —> bb/i+yU~ are displayed in Fig. [6] for the generated 
(hatched) and reconstructed (full color) events without any cut applied . 

The preselection requirement on muon pair (cut 1, Tabled]) reduces approximately to same 
fraction of signal (~ 62%) and background (Z~ 59%, tt ~ 67%) events. The requirement on 
jet pair (cut 2) leaves almost untouched tt (~ 63%) reducing Z to ~ 54%, because of the more 
energetic -Prjet spectra of the top decay [7]. Unfortunately for the same reason this cut reduces 
the signal sample significantly to ~ 30% , see Table [71 

The mean jet transverse momentum < -Prjet > is ~ 26 GeV for the h boson, and ~ 36 
GeV for Z decays, see Fig. [7] [7]. Then, as stated in Sec. 13. 3[ the requirement that two b-jets 
be identified will not enhance the signal over this specific background. This requirement will 
suppress more the signal than the Z background, since the two b-jets of this background events 
are usually more energetic than those accompanying the signal. According to the study in Sec. 
16. 2^ only one jet is required to be tagged as a b-jet. 

After the three preselection cuts (1-2-3) the signal sample of h or A is reduced to ~ 10%. 
The backgrounds are reduced to ~ 17% for the Z, ~ 53% for tt, and ~ 28% for ZZ. 

The tt veto selection is designed to exploit features of t decays. The tt sample is char- 
acterized by a transverse missing energy E^f^^^, which, due to neutrinos, is larger than in the 
h/A sample, as shown in Fig. [HI The E'™^^ cut (cut 4, Table [5]) is extremely effective, reducing 
tt to ~ 14 % of the original sample while keeping the signal at ~ 9 % almost unchanged and 
the Z background at ~ 14 %. The two cuts (5-6), related to Pt^i distributions, reduce the tt 
original sample to ~ 8 %, but are little effective on the Z sample (to ~ 10%) and the signal 
samples of h and A (to ~ 8 %). A further reduction of the tt background to 1.4 % is obtained 
by applying cuts (7-8) on Prjet- At this stage, the Z sample is reduced to ~ 5.8 % while the h 
signal as well as the A signal is reduced to 6.0 %, see Fig. [HI 

The effect of cuts (1-8) is summarized in Fig. [TU] which shows the distributions of the 
reconstructed invariant mass for signal (h and A) and (weighted) background (Z, tt and 

ZZ added up) events. The h/A signal (light blue) is clearly visible on top of the remaining 
background events (Z, tt and ZZ added up, dark brown). 

To evaluate the significance of observation of a signal, the mass window (Eq. \^ is then 
selected, final selection (cut 9, Table [3]). At the reference point {nix = 110.31 GeV, Fa = 4.28 
GeV and mh = 110.00 GeV, Fh = 4.20 GeV), Eq. [5]implies to look for the A signal in a window, 
^A, of 3.394 GeV around a corrected mass value, mA^, of 109.490 GeV, and for the h signal in 
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Figure 6: Without cut. Distributions of the number of muons iV^ (top) and of the transverse momentum Pt/j 
(bottom) for generated (hatched) and reconstructed (full color) events are plotted without cuts applied: 
a) bbh — + bb/i+/i^ events at the reference point (tan/? = 45, Wh = 110.00 GeV)(left, dark blue); b) bbZ 
bb/i+/i^ events (right, yellow). All distributions are normalized at j Cdt = 300 fb~^. Entries are per bin width 
of 1 (top plots) and of lO'^ MeV (bottom plots). 




PijeJMeV] 




Figure 7: After 1-2 cuts. Distributions of the reconstructed jet transverse momentum Pxjot are plotted after 
cuts 1 and 2_ (Table E]) for : 

a) bbh hhfi^^" events at the reference point (tan/3 = 45, nih = 110.00 GeV), (left, dark blue); b) bbZ 
bb/x^/i" events (right, yellow). All distributions are normalized at j Cdt = 300 fb~^. Entries are per bin width 
of 103 MeV. 
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Figure 8: After 1-2-3 cuts. Distributions of the reconstructed transverse missing energy E^^^'^ are plotted 
after cuts 1, 2 and 3 (Table E]) for : 

a) bbh bb/^+/i^ events at the reference point (tan/3 — 45, mh = 110.00 GeV), (left, dark blue); b) tt 
hh^'^ IJL~vv events (right, gray). All distributions are normalized at j L At — 300 fb^^. Entries are per bin width 
of 10^ MeV. 



a window, of 3.369 GeV around m^°" = 109.180 GeV (see [7] for the others values of tan/3, 
TTiA and mh). As for the probabihty requirement f = 2, the windows and are increased 
by a factor 2 to obtain the effective selection window used in Table O 

This selection reduces drastically the background samples, the Z to less than 0.13%, the 
tt to 0.17 % and the ZZ to a negligible 0.24 %. The event signal surviving in both A and h 
channels is ~ 4.2 %. 

The last reduction (cut 10) on muon isolation does not change significantly the previous 
results. Indeed this cut is expected to be effective in the data against the heavy flavour QCD 
background, and cannot be tested properly in Monte Carlo. 



Process 


]\^exp3oo 




Nlcut 


N2cut 


Nscut 


bbh hhfi^fi^ 


73500 


76517 


47099 


22969 


7907 


bbA — > hh^'^^~ 


72900 


74965 


46247 


22693 


7919 


bbH — s> hhfi~^fi~ 


486 


600 


394 


222 


95 


bbZ — > bb/i+/i~ 


6836700 


3314000 


1945387 


1178864 


579751 


tt — i> hhfi~^fi~h'u 


1713420 


1806437 


1217706 


1137780 


953928 


ZZ bb/i+/i~ 


33819 


97244 


46487 


38600 


27744 



Table 7: Preselection. Signal/Background process, number of events expected at JCdt ~ 300 fb^^, N^'^p^"", 
number of Monte Carlo events generated, N'^^, and after three preselection steps, Nicut, N2cut and Nscut (Table 
[5]). The signal is evaluated at the reference point (tan/3=45 , tua — 110.31 GeV, mh — 110.00 GeV). Cuts 1 to 
3 are described in Table [5] 

We calculate that the significance at JCdt = 300 fb~^ is ~ 56, scaling down to ~ 18 
at fCdt = 30 fb~^ and 10 at JCdt = 10 fb~^. As discussed in Sec. 13. 3[ the process of Z 
production accompanied by two light jets, which has not been simulated, can be an additional 
significant background. Its importance has been estimated from the relative cross-section and 
fake b-tagging probability. When this background is taken into account, we have estimated 
that the significance is lowered by ^ 25%. 
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Figure 9: After 1-8 cuts. Distributions of the reconstructed invariant mass M^"^^ are plotted after cuts 

1-8 (Tablep for: 

a) bbh hh^~^^~ events at the reference point (tan/3 = 45, mh — 110.00 GeV) (left, dark blue); b) bbZ 
bb/x^/i" events (right, yellow). All distributions are normalized at JCdt = 300 fb^^. Entries are per bin width 
of 1.5 • 10^ MeV. 



Process 


Nscut 


N4cut 


Nscut 


Ngcut 


Nrcut 


Nscut 


bbh — s> bbyU+yU^ 


7907 


7099 


6234 


5575 


5085 


4636 


bbA bb/i'^/i~ 


7919 


7096 


6183 


5575 


5118 


4595 


bbH bb/i"^/i^ 


95 


78 


66 


54 


50 


44 


bbZ —>■ bb/i+yU^ 


579751 


478106 


397030 


334609 


226414 


193334 


tt —>■ hhfi^ ^^uu 


953928 


256662 


190672 


143746 


41954 


25567 


TIL bb/i+yU~ 


27744 


24516 


19102 


15957 


10926 


8574 



Table 8: tt veto selection. Signal/Background process, number of Monte Carlo generated events after the 
preselection cut 3 and cuts 4 to 8. The signal is evaluated at the reference point (tan/3=45 , toa = 110.31 GeV, 
TOh 110.00 GeV). Cuts 3 to 8 are described in Tabled 

We conclude that if = 110.00 GeV and consequently mj^ = 110.31 GeV there is a high 
probability for these bosons to be discovered at the beginning of data taking. 



8 Search in the MSSM plane 
8.1 bbA hhfi+fi-, bbh bbfi+fi' 



A search for neutral Higgs bosons h/A has been performed in the (mA, tan/3) plane, varying 
tan/3 between 15 and 50, in steps of 5, and mh between 95 GeV and 125 GeV, in steps of 2.5 
GeV. For each point the statistics reached leads to weight factors for h and A boson close to 
unity. 

The analysis described in Sec. [7] has been repeated for all 104 simulated points and the 
results are reported in Ref. [7]. 

The use of an asymmetric selection window, centered at m^p^ (see Sec. 17.11) . has been tested 
for masses below 100 GeV, with a view to exclude the Z events, however without obtaining a 
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JT 1 ULcoo 


Scut 


9cut 


lOcut 


UUli ^ UUjlX jlX 






3165 






fj X tJ X 


31 1 


hhH hh;; + ;;" 

UU±± ' \J VJ Lil jLfc 


44 






hhZ — > hh;; + ;;" 

UUZJ ' yJULij Li 




44Qfi 


4206 


ft, — llVl//^// 777/ 


25567 


3148 


2969 


TIL bb/i+/i^ 


8574 


260 


234 



Table 9: Final selection. Signal/Background process, number of Monte Carlo generated events, after cut 8 and 
successive cuts 9 and 10. The signal is evaluated at the reference point (tan/3 — 45, tua — 110.31 GeV, mh = 
110.00 GeV). Cuts 8 to 10 are described in Tabled] (see text for cut 9 window). 

significantly improved result (~ 10%). A symmetric window as described in Sec. 17.11 was thus 
applied to all masses. 

The search significance for the h/A neutral boson is shown as a function of rriA up to highest 
allowed value of mh, in Fig. [H] for all scanned values of tan/5 and three luminosity, JCdt = 300 
fb~^ ('5h°A), 30 fb-i (^^%) and 10 fb'^ (^^%) [7]. The values for the two lower luminosities 
were derived from the first one, which corresponds to the highest statistics.. 
One should note that large h/A masses are penalized by a small cross section, thus implying 
a lower significance, while the masses near to mz suffer from the difficulty in disentangling the 
neutral Higgs boson signal from the Z background. 

The best mass range for an early discovery of h is between 100 and 120 GeV at any given 
tan/?. If tan/3 > 30 a large range of masses is accessible to discovery even after the first year 
of data taking. More integrated luminosity, between ^ 30 and 50 fb~^, is needed for tan/? 
between 30 and 20. The discovery at tan/3 = 15 demands a luminosity of ~ 150 fb~^, making 
the exploration of this region possible only after a few years of data taking. 

These considerations are summarized in Fig. [12], where the minimum integrated luminosity 
/ Cdt, demanded for a 5 a discovery of the h/A neutral Higgs boson, is plotted as a function 
of mA up to highest allowed value of mh, at given tan/3 values. With a /£dt ~ 10 fb~^, 
corresponding to one year of data taking, most of the masses are accessible if tan/3 > 30. More 
integrated luminosity is needed for tan/3 = 20 and tan/3 = 15. Low masses need as well more 
luminosity in order to extract the evidence of a signal from the most copious Z background. 

Discovery contours in the (tan/3, plane are shown, in Fig. [121 in different j Cdt scenarios 
for a significance of 5 (discovery, on the left) and of 3 (on the right). The latter can interpreted 
as the contour region for an early indication of a signal or, in case of negative search, for its 
exclusion. 

8.2 bbH hhfi+fj,- 

To increase the discovery potential for a neutral Higgs boson discovery in the region up 139 
GeV a search for the neutral Higgs boson H has been performed in the (tan/3, mA) plane inside 
the limits used for the h/A search (Sec. 18.11) . This search required a separate analysis due to 
the higher mass region involved. 

Due to the extremely low cross section of the H in this mass region, for a few masses with 
a cross section times branching ratio crH-Br^+^- smaller than 0.01 pb, common values 
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Figure 10: Distributions of the reconstructed fJ-^ IJ-^ invariant mass, M™^ , for signal and backgrounds events, 
after the selection cuts 1-8 (Table [5|) at the reference point (tan/3 = 45, ttt-a = 110.31 GeV, mh = 110.00 
GeV). The two distributions are normalized at ^Cdt— 300 fb^^. The h/A signal (light blue) emerge over the 
background (Z, tt and ZZ) (dark brown). Entries are per bin width of 1.5 • 10^ MeV. 

of center and selection window were chosen, corresponding to the mean values in that mass 
region. This procedure doesn't affect any result due the narrow spread of values for the mass 
mn and the natural width Fh, implying that the window value is essentially dominated by the 
experimental resolution. 

The discovery of the boson H demands a high integrated luminosity, j Cdt ^ 300 fb^^, and 
it would be possible at masses corresponding to mh = 122.50 GeV and 125 GeV for values of 
tan/3 > 30. At tan/? = 25, only a narrow range of mn masses around 134 GeV (corresponding 
to mh = 125 GeV) is accessible. Lower values of tan/3 are excluded even at high masses. 

In conclusion, a discovery of H boson is possible only in a few points of the parameter space, 
for a value of rriA around the maximum value of mh, at the ultimate luminosity expected at 
the LHC. 

8.3 Combined search for bbA — > bb/x+//~, 
bbh — > bb/i,+/Lt~, bbH — > bb/x+/x~ 

The results of Sec. 18.21 on the H boson search were combined with the results from the h/A 
search (Sec. 18.11) . To this purpose, the H analysis was repeated with a window not overlapping 
with the h/A window, thus avoiding the double counting of background events. This work was 
performed only for high mh values, where the cross section times the branching ratio 

(Sec. EM has values > 0.01 pb [7]. 

In Table ITO] the significance of the exclusive search for the h/A boson (5*;^°^' a) 
and the corresponding significance including the H search H' a h a h) 

/£dt = 300, 30, 10 fb-^are shown. 
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Figure 11: Search significance for a h/A neutral Higgs boson, as a function of ttia up to the largest allowed 

value of mil in three different data taking scenarios, j Cdt = 300, 30 and 10 fb~^ (S is the number of h/A signal 
events, B is the number of background events). On the left the results for tan/3 = 50, 40, 30, 20, and on the 
right the results for tan/3 = 45, 35, 25, 15. The data are listed in Ref. [7]. 



22 



Lum {5a) vs m 
• tanp = 50 
T tanp = 40 
tanp = 30 
■ tanp = 20 



%0 95 100 105 "110 "115 '120 125 130 135 140 

[GeV] 



-"300 

J25O 



200 
150 
100 
50 



Lum {5a) vs m 
• tanp = 45 
T tanp = 35 
A tanp = 25 
■ tanp = 1 5 



f i a i » i i l l Bi i # i #»T i I * 



%0 95 100 105 '1f0'1l5 "120 125 130 135 140 

rriA [GeV] 



Figure 12: Minimum integrated luminosity JCdt demanded for a 5 cr discovery of the h/A neutral Higgs boson 
as a function of ttt-a up to largest allowed value of mh- On the left the results for tan/3 = 50, 40, 30, 20, and 
on the right the results for tan/3 = 45, 35, 25, 15. The blue horizontal line is JCdt = 10 The data are 

listed in Ref. [7] . 
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Figure 13: 

Discovery potential for a neutral Higgs boson h/A of mass mA decaying to /x^/i^, accompanied by two b-jets, 
in the mh — max scenario (Sec. [2]), as a function of m^: contours are drawn for a search significance = 5 

(left) and — 3 (right), with an integrated luminosity of JCdt = 300 (top), 30 (center) and 10 (bottom) 

fb-i. 
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The additional H search contributes to the MSSM Higgs sector discovery for tan/3 > 15 at 
the ultimate luminosity. Otherwise the contribution of this search is negligible, for a value of 
mA around the maximum value of rrih, as expected from the low cross section, and the large 
number of background events (being mn close to mz). 
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Table 10: Significance of A/h and A/h/H searches. The significance values for the two searches, ^h^A' "^h^A' 

•^h" A^^'i '5'j^°A H' H' "^h^A H'^* /£ di = 300, 30, 10 fb~^ , respectively, arc givcn for a Set of tan/3 values 
and two values of the mass of the lightest neutral Higgs boson, to{J°™. The mass nin of the corresponding H 
Higgs boson is also noted. 



9 Conclusions 

The possibility of the discovery of the MSSM h/A bosons in the region of high tan/3 ( larger 
than 15) and mass close to 100 GeV has been investigated by exploiting the decay of the neutral 
h/A boson into two muons, h— > //"'"/i" and A— > n'^/i' , accompanied by two b-jets. This region 
is also accessible by charged MSSM Higgs boson H='= decays. 

For this purpose an analysis using full detector simulation has been performed. Monte Carlo 
events have been generated for a center-of-mass energy ^/s= 14 TeV through the ATHENA 
interface (v.9.0.4), while the ATLAS detector response has been simulated using the GEANT 
program through the ATHENA interface (v.10.0.1). 

The results described in this paper show a well defined possibility for the discovery of a 
neutral Higgs boson in a region traditionally difficult due to the presence of the Z resonance. 

This is achieved thanks to the high resolution performance of the ATLAS detector, namely 
of the muon spectrometer and the inner detector, together with the high b-tagging capability. 
For completeness the search of H— > has been explored in the same mass region. 
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The discovery of a neutral MSSM boson looks possible in a mass range of 100 to 120 GeV 
at tan/5 > 15, with an integrated luminosity j Cdt = 10 fb~^, which corresponds to one year 
of data taking. A warning should be made however because of unaccounted uncertainties in 
background and other systematic errors that cannot be evaluated precisely at this time. 

With a view to perform this analysis on real data a method to subtract the main contributing 
background of Z boson decays Z — has been suggested by us [8]. This method mainly 
relies on experimental data with limited Monte Carlo corrections. The procedure is based on 
the use of a control sample of Z bosons decaying to electrons, and does not depend on complex 
theoretical calculations nor on their implementation in Monte Carlo. 
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A Appendix: 

subtraction of the bbZ — > bb/x+/x~ background 

Following the method outlined in Sec. 16.31 and based on the universality of the lepton coupling, 
the number of background events bbZ bb/i"'"/i~, Nz^^+^-, is given by the number, Nz^e+e-, 
of bbZ— >e"'"e~ events collected in the same real data sample. In both cases the number of events 
is the one counted within a window of the two-lepton invariant mass distribution. The method 
has been fully proved with the simulation described in [8]. 

The ratio Ni-^^+ ^- / Nz^^+e- has been shown to be a regular function of the dilepton in- 
variant mass in the region of interest between 97.5 GeV and 140 GeV once possible sources of 
differences between e"'"e~ and data have been taken into account, and to be close to one. 

We like to recall here these differences and their impact on the above ratio. 

Clearly, depending on whether the final state contains electrons or muons, one has to provide 
for the differences in the detector response that is for the acceptance and resolution of the 
electromagnetic calorimeter and the muon spectrometer. 

Since our samples are bbZ fi~^fi~ or bbZ — >■ e"'"e~, in those events where a semileptonic 
b decay occurs with the same flavour as one of the Z decay products, more than one invariant 
mass combination is possible, namely three (only opposite charge lepton pairs are considered 
for invariant mass combinations). Often, the value of the invariant mass associated to this 
combination of one lepton originating from Z and another one from b decay sits far from 
mz (and thus has no effect in the interested region). However the higher efficiency in muon 
detection increases the number of fake combinations in the muon sample. 

Less obvious are the corrections required by the different inner bremsstrahlung in the two 
samples. The inner bremsstrahlung (IB) is the emission of photons, 7brem; near the Ze"'"e~ or 
Z/i"*"/!" vertex. The presence of such photons changes the kinematic configuration of the decay, 
in particular the lepton 4-momenta. The impact of this effect on the Higgs mass resolution with 
the ATLAS detector has been studied in Ref. [33]. Here we discuss its impact on the number 
of detected events as a function of the dilepton reconstructed invariant mass in the region of 
our search. 
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At the generator level in the electron (muon) sample 85.2% (91.6%) of the events are without 
IB photons, 13.7% (8%) have 1 7brem, 0.94% (0.2%) have 2 7brem- The corresponding average 
transverse momentum of the radiated photons, without any cut except those mentioned in Sec. 
m is close in the two samples, < P^^""^ >= 16.18 GeV for the e"'"e~ sample and < Pr^'"'™' >= 
16.43 GeV for the sample. These IB photons are mainly contained in a cone with an 

opening AR = 0.15 around a muon track and AR = 0.25 around an electron track. By 
summing up the IB photon 4-momentum to the close lepton 4-momentum one obtains that the 
ratio Nz^^+fj,- /Nz^f,+c- stays the same within ±12% for dilepton invariant masses between 100 
and 140 GeV. It would otherwise show a spread of ~ ±30%. 

At the detector level muons, electrons and photons involve different detectors (muon spec- 
trometer and electromagnetic calorimeter) implying different momentum resolution, angular 
acceptance and efficiency. In this analysis we used identification criteria as follows. 

A muon is reconstructed as a track, in both the inner detector and the muon spectrometer 
(so named combined reconstruction), and enters the analysis if Pt > 10 GeV and \rj\ < 2.5. 
The average momentum reconstructed is < Pp >= 37.83 GeV. The reconstruction efficiency 
is higher for muons than for electrons. 

An electron is reconstructed as a cluster in the electromagnetic calorimeter associated with 
a track inside Arj = 0.025 and A0 = 0.05, and enters the analysis if Pp > 10 GeV and 
|?7| < 2.5. The best performance is achieved when the energy is measured in the electromagnetic 
calorimeter (with a cluster size of 5 x 5 cells) and the angles {6, 0) in the tracker. The average 
momentum reconstructed is < Pp >= 42.61 GeV. 

A photon is reconstructed as a cluster in the electromagnetic calorimeter not associated 
to any track inside Arj = 0.025 and A0 = 0.05, and enters the analysis if Et > 10 GeV and 
\rj\ < 2.5. A cluster size of 5x5 cells is used to measure the energy. The value chosen for the 

threshold ensures a good photon reconstruction efficiency. 

In this analysis however a photon's 4-momentum is summed with that of the lepton if 
the photon is reconstructed inside a cone AR < 0.15 around an electron (e^e~ sample) or 
AR < 0.25 around a muon {fi~^fi~ sample). As a result, although the number of photons of 
any origin is more important in the e+e~ than in the sample, the number of photons 

reconstructed separately is equal (~ 1%) in the two samples. The effect of summing the photon 
and lepton 4-momenta shows up in the distribution of the ratio N^+^- /N^+e- as a function of 
the dilepton invariant mass when only events containing an IB photon at the generator level 
are considered. It consists in a shift from the low-mass tail towards the Z mass. 

Based on the whole simulation we could conclude [8] that the ratio N^+f^- /N^+e- is ~ 1.2 
and does not depend, in the region of our search and within the simulation statistics, on the 
two-lepton invariant mass. 
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